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Abstract 



Detailed study of multifractal characteristics of the financial time series of 
asset values and of its returns is performed using a collection of the high fre- 
quency Deutsche Aktienindex data. The tail index (a), the Renyi exponents 
based on the box counting algorithm for the graph (d q ) and the generalized 
Hurst exponents {H q ) are computed in parallel for short and daily return 
times. The results indicate a more complicated nature of the stock market 
dynamics than just consistent multifractal. 

I. INTRODUCTION 

Comprising the financial time series into a unique description on different time scales 
by means of the multifractal formalism is an attractive perspective both for fundamental as 
well as for practical reasons JT|. As a constructive, though also critical, contribution towards 
this aim in the present paper we systematically review the related basic characteristic of the 
financial time series using the large collection of the high-frequency Deutsche Aktienindex 
(DAX) data with over 10 6 data records within the time period 28. Nov 1997 4 30. Dec 1999 
0, as well as the daily data since 29. Sep 1959 4 4. Sep 2001 J| (over 10 4 data points). We 
thus estimate the tail index for the distribution of DAX returns, the Renyi exponents for the 
DAX index and return's graphs, as well as the corresponding generalized Hurst exponents. It 
should be stressed that contrary to all the previous related studies the fractal exponents for 
the graphs are here computed directly, using the (2-dimensional) box counting algorithm, 
independently of the Hurst exponent computations. 

Let S(t) represent an index value as a function of time. The rate of return within the 
time lag r (the return time) is usually denned as 

S(t + r)-S(t) 

GAt) ~ W) ' {) 
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The logarithmic returns 



G T (t) = \nS(t + r) - \nS(t) = In (1 + G' T {t)) ~ G;(t) (2) 

are more often used when analyzing the financial time series. They are additive with respect 
to the return time, 

G T1 (t) + G T2 (t + r l ) = G T1+T2 (t) . (3) 

This implies translational symmetry of returns and scaling invariance of the asset price 
(return is identical for prices in any currency). The logarithmic returns are indistinguishable 
from normal returns if <C 1. In general, the normal returns ([I]) are slightly bigger in 

absolute value which leads to fatter tails due to the inequality x > m(l + x) for x G (0, 1). 
Simple index differences have also been investigated but they are dimensional quantities [|J . 
Throughout this paper by the index x t at time t we will understand the logarithm of its 
value, x t = \a St- Hence, the returns G T are increments of the index x t , just as the white 
noise (WN) random variables are increments of the Brownian random walk (BRW). The 
index x t can be viewed as integrated returns 



E 

k=ti 



x (t 2 +i)r = x tlT + G T (kr) , (4) 



where xt lT is the integration constant. To compare various returns the normalized returns 
are introduced that have zero average and unit variance 

G T (t) - (G T (t)) t 
gAt) = a(G T (t)) ' (5) 

where (. . .) t denotes averaging over time variable and a(. . .) is the standard deviation. 

In the following section the probability density function (pdf) of the DAX returns is 
discussed in order to demonstrate that it shares essentially all the global characteristics 
identified so far in the literature for other indices, though some small quantitative differences 
do of course take place. Section III is devoted to study self similarity of the DAX index 
and of its returns using the (2-dimensional) box counting algorithm. The generalized Hurst 
exponents are analyzed in Sec. IV. Summary and conclusions are given in the last section. 



II. DISTRIBUTION OF RETURNS 

The basic property of normalized returns that has been intensively investigated in recent 
years is their pdf (obtained by constructing histograms from the empirical data). The old as- 
sumption that returns follow the uncorrelated BRW with independent identically distributed 
(iid) Gaussian random variables |J has been challenged by Mandelbrot and by Fama ||[7|. 
Mandelbrot analyzed a relatively short (about 2000 data points) time series of cotton prices 
and he reached conclusion that returns have Levy stable symmetric distribution with Pareto 
(power like) fat tails 

This result has been refined in recent years by the Boston group |10| , pT|| using over 10 6 
high frequency data points of the S&P500 index. They have found clear departures from the 
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Levy stable distribution. First, suggesting the exponential truncation of fat tails for returns 
> 6er, where a is the standard deviation. Later, they found a power law fit with the tail 



index exceeding the Levy distribution [fTI |. The last result was confirmed by other groups 
for distributions of returns of other indices, including the FOREX data |12| . 

In this Section we investigate histograms of the DAX returns for high and mid-frequency 
data. In Fig. ^ we plot the minutely data histogram (r = 1 min.) for normalized returns 
(0). The normalized Gaussian distribution (dotted line) as well as the best fit for the Levy 
stable distribution (dashed line) are also plotted. It is clear that the high frequency returns 
are neither Gaussian nor Levy stable. The best fit for Levy distribution within the central 
region g G [—gi] +gi], with g I ~ 5 gives for the Levy parameters a = 1.16 and 7 = 0.286. It 
is interesting to notice that this fit is very stable with respect to the range of the variable g: 
one gets identical results within about 1% for g>j = 1, 3 and 10. The tail index (a) obtained 
in this way is too small to fit well for g > 2. Also, it is smaller than the value 1.4 for S&P500 
reported in |TT . 



Here, one should notice one exception: the central bin contains abundant number of 
data points (and it was not taken into account for the fits). This phenomenon we call 
the zero return enhancement and it is typical for high frequency data, especially for single 
stocks, indices that do not contain too many stocks and for not too large stock markets. 
For short time lags it happens that there are no transactions, the index value is constant 
and the return distribution is enhanced for the single return value g = 0. For 15 sec. 
DAX returns this effects gives about factor 2 enhancement for the central bin and still 
considerable enhancement for r = 1 min. However, for much larger S&P500 index (and 
larger NY stock market) this effect is hardly visible even for minutely data. Because of this 
effect and discreteness of returns, we prefer the least square fit adopted for Levy distribution 
I3| within the central region of the histogram, with the zero return bin neglected, instead 



of using the scaling form of the return probability p(0) (see also |10| and Sec. IVA in JTT 



The log-log plot corresponding to Fig. [I] is given in Fig. |5|(A). For comparison the dashed 
line with slope corresponding to the tail index a = 2.4 is also plotted. We can see that both, 
negative and positive returns give approximately the same tail index close to the value given 
above. This value considerably exceeds the range of the Levy stable distributions. Our 
estimates give a ~ 2.5, in particular a = 2.4 for positive and a = 2.6 for negative returns. 
Here, one can observe an effect that is reverse to what was found in other markets: the 
negative return tail has over 8% bigger exponent than the positive tail. Quality of the linear 
fit for the n data points with coordinates (xj, yi) is usually measured by the linear correlation 
coefficient r (also known as the Pearson's r), where r e (— 1,+1). For r = 1 the fit is of 
course perfect. Here, in both cases the correlation coefficient r > 0.996 within the range 
2 < g < 20. The resulting tail index, a, is considerably below the S&P500 result (a ~ 2.8), 
reported in [|ll|] (though, in both cases well above the Levy stable limit a = 2). In practical 
terms the DAX develops somewhat fatter tails which reflects its more sizable fluctuations 
(more "rare events"). 

For g > 20 (or log 2 g > 4.5) the bins of the histogram contain very few data points (0-4-3 
for reasonable bin width that is ~ 0.1 for Fig. 0(A)) and any regression for those rare events 
is rather doubtful. Hence, any fit for the very far tail (g > 20) truncation seems to be not 
well justified even for the largest sets of the high frequency data available today. The same 
results have been also obtained for r = 15 sec. In summary, the high frequency DAX return 
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distribution appears to be quite complicated consisting of the zero return enhancement in 
bin with g = 0, Levy like distribution in the central region (\g\ < 2) with a ~ 1.3, inverse 
power like behavior with a ~ 2.4 for the tail (2 < g < 20) and not easily quantifable 
behavior for very large events (g > 20). 

The above picture, with two different a's for the central region and the tail, seems to 
be consistent with recent suggestions that the return distribution is a union of different 
distributions for returns in "normal" day time and "rush hours", like soon after opening, 
before closing and at 14:30 in Europe |l4]-^7[]. In this paper we however do not extract special 
subsets of our financial time series that would screen out such effects as they constitute an 
essential element of the stock market dynamics. Doing so would also shrink the data set 
which even as a whole seems to be of relatively modest size (see e.g. flTj||). Furthermore, 
in addition to trends in the special subsets mentioned above there exist some other trends 
on various time-scales which are not connected with the calendar time but are associated 
with the internal market dynamics (bubbles and crashes, see |iTS| , |2Ti|| for a recent discussion). 
Hence, throughout this paper we analyze the complete DAX data set available. 

Increase in the return time, r, leads to the increase of the parameter a of the fitted Levy 
curve. In FigsfJH and Figs.|5](B,C) the same plots as in Figs.[l]j5|(A) are displayed but for 
return times r = 10 min. and r = 1 hour, respectively. In both cases we have good fit of 
the Levy distribution in the central region, g G (— gj, +gi), gi — 3. The Levy parameters 
a and 7 slowly grow from the values a ~ 1.2 and 7 ~ 0.29 for the minutely returns up to 
a ~ 1.4 and 7 ~ 0.39 for the hourly returns (see Figs. 0,0,0). At the same time the tail 
index outside the central region grows from a = 2.4 for the minutely returns up to a ~ 3.5 
for hourly returns (see Figs. 0(A-C)). 

For return times longer than r = 60 min. the number of data points is decreasing below 
10 3 which is not sufficient to precisely estimate the tail behavior. Hence, for mid-frequency 
data we use much longer data set with historical daily data to have reasonably long time 
series (> 10 3 ). The return time here, r ~ 500 min., is over 3 orders of magnitude longer 
than for the shortest available (r = 15 sec.) data. Again, the central region where the 
histogram can be well approximated by a Levy stable distribution is: \g\ < gj ~ 3. The 
best fit for the Levy parameters in this region (also, quite stable with respect of the changes 
in gj) is: a ~ 1.7 and 7 ~ 0.385 (see Fig. f|). In the lower panel, we have added 40-fold 
magnified tail to show that even for the largest available set of data the tail behavior cannot 
be reasonably estimated already for g > 5. For longer return times the situation is even 
worse. For r = 30 min. the tail index seems to be stabilized at the hourly return value 
a ~ 3.5 (Fig. 0(D)). In all cases it is greater than the maximum value allowed for the Levy 
stable distributions (a < 2.0), in agreement with what was found for other indices. 

Good fit of the Levy stable distribution within the central region is consistent with the 
old result obtained by Mandelbrot for a data set of quite modest size However, change of 
the Levy parameters with the return time (r), as well as differences in the value of the Levy 
parameter a for the central region and for the tail in each case clearly shows the complexity 
of the return distributions and their incompatibility with a plain Levy curve. Also, from 
lower panel of Fig. (| one can see apparent difference in the tail behavior of the normalized 
return distributions and the normalized Gaussian distribution. 

The behavior of rare events (i.e. for tail with normalized returns g > 20) cannot be 
reasonably estimated for return times r > 1 min., at least using the standard histogram 
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building and fit method, due to the insufficient number of such events in the present day 



financial data available. Similarly as for the S&P500 fil l, even though on short time scales 



the DAX return distributions are clearly not Levy stable, in the range from r = 1/4 min. 
up to r = 1 day we do not observe convergence of the pdf to the Gaussian (see Figs. []]-§]). 

For the standard (Lindeberg) version of the Central Limit Theorem (CLT) the random 
variables must be iid and the stochastic process must be stationary. This condition is not 
satisfied for financial time series where, in particular, various short as well as long range 
cycles are present. In particular, empirical data display clear nonstationarity as can be seen 
from the pattern of the recurrence plot that is given for 10 min. returns data series in Fig. 
Here, the relation between closeness in time and in phase space is displayed using the Takens 
phase space reconstruction method. The embedding dimension is taken E = 5, but its exact 
value is not important in this the recurrence plot is not sensitive with respect to 

this parameter. The normalized neighborhood size was set e = 0.2, but similar structures 
were obtained for smaller e (see e.g. |2T| Sec. Ill and [^,^] for more details), as well as 
for the embedding dimension up to E = 9. Recently, Bouchaud et al have shown that 
for a special model with correlations similar as in financial series a slow convergence to the 



Gaussian pdf can be proven and it is visible for return times r > 1 day (see also |TT 



III. SELF SIMILARITY 

Distribution of returns gives us only a small part of the whole information included in the 
time series. In particular, the time order is completely neglected. To account for the time 
dependence one should analyze the whole (2-dimensional) graph of the functions x t = x(t) 
and g T (t) (D G in the notation of 0). Although, the projections of the data points at the 
time and at the value axis are pseudofractal sets with zero Hausdorff-Besicovitch dimension 



and their pseudofractal scaling exponents are equal to one [25], the graph in principle can be 



a regular fractal as is the case for the WN and BRW. The fractal properties of returns have 



been suggested long ago by Mandelbrot pq| . More detailed analysis has been published by 
Evertsz et al [I2,27[|. However, it should be stressed that their analysis was based on the 
calculation of the generalized Hurst exponents (see Sec. IV). 

To investigate fractal properties of DAX index and its returns we apply directly the 
standard box counting algorithm for computation of the Renyi box counting exponents 
defined by [^8 



i h b5ffl 3h m, (6 ) 

q 1-q N-.OC In N a^oo In A w 

where N is the total number of "boxes" (bins), Pi is the part of the "mass" (i.e. fraction of 
all points) contained in the i-ih box. 

Computation is performed for both graphs, x(t) and g T (t). To the best of our knowledge 
this is the first direct calculation of the fractal dimensions of the graphs for any financial 
index and its returns. Usually, the self similar properties were inferred from the Hurst 
exponents or from the rescale range analysis, see e.g. [^,^]. As will be shown in the 
following sections these two methods give different results, especially for the returns' graphs. 
In Fig. [7] results of the calculations of d q exponents are displayed for the logarithm of DAX 
(xt). The calculations for time lags r = 1,10,60 min., 1 day and for q = 0.5,1,2 and 4 
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are shown. In fact, the same results were obtained for the plain index (St). This supports 
the view that the graph is a real fractal, not just a pseudofractal as the exponents for 
fractals are invariant with respect to the homeomorphic transformations. In all cases one 
gets quite reasonable linear fits. The correlation parameter (r), as can be expected, is 
lowest for r = 60 min, where the number of original data points is the smallest (below 10 4 ). 
However, even in this case we have r > 0.998 (for all other cases we have r > 0.9994), 
quite a reasonable value. The linear scaling extends through 6 -j- 9 points (binary orders 
of magnitude) at the log-log plot which is also a standard range for numerical estimates of 
fractal dimensions. The resulting Renyi box counting exponents are given in Table |. Even 
though, the range of the return time (r) and the range of the index q are relatively large, 
all values for the exponents d q are within the interval d q € [1.3, 1.4], well below the value 
characteristic for the BRW (in all Tables the last digit is not the significant digit). This 
suggests that we have a self similar and close to monofractal curve for the DAX index. In 
addition, one can observe a general tendency that the index d q is slowly decreasing for longer 
time lags except for r = 1 day, where another set of data was used. 

Similar analysis applied to the normalized DAX returns leads to the result displayed 
in Fig. H Here, the quality of the linear fit in the log-log plot is not as convincing as in 
the previous case (Fig. |7|) of the index itself, especially for the smaller g-values. Extracting 
nevertheless the corresponding linear best fit (r > 0.994) coefficients results in the values 
for d g which are more dispersed and spread within the range d q = 1.4 1.7, depending on 
q as well as on the return time, r. 

Our results thus indicate that, within a few percent, which is comparable to the es- 
timated numerical accuracy, the g-dependence (and r-dependence) of d q is rather weak. 
Hence, the corresponding graphs can be viewed as close to monofractal and this indication 
is especially suggestive for the index whose scaling properties are quite convincing. It is 
also worth to notice that the corresponding d q values for the index are systematically be- 
low the BRW value (1.5). This may reflect an important element specific to the financial 
dynamics. Somewhat less transparent is the related scaling behavior of the returns, though 
this behavior also carries some information about the underlying dynamics and is therefore 
shown here. Assuming however that some scaling approximately applies in this case as well, 
it seems natural that the above estimated bounds limiting variation of d q for returns point 
to values which are systematically bigger than for the index, as the former is the differenti- 
ated index (see eq.(P). While this difference equals 0.5 for WN and BRW, in our case the 
difference is on average significantly smaller. To avoid numerical artifacts we have repeated 
all calculations with logarithm of the index and we have obtained the same results within 
the numerical accuracy. This shows that the results are invariant with respect to nonlinear 
transformations of the data. As the fractal dimensions and Renyi box counting exponents 
are related to the Hurst exponents we will return to this issue in the subsequent section. 



IV. GENERALIZED HURST EXPONENTS 



The Hurst exponent (Hi) has been introduced long ago in the "rescale range theory" 
for measurements of the Nile flooding and drought amplitudes pi |. The generalized Hurst 
exponents, H q = H(q), for a time series g(t) (t = 1,2,...) are defined by the scaling 
properties of its structure functions S q (r) PJ 
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S 9 {T) = (\g(t + T)-g{t)\*) 1 J q ~'r JI ® (?) 

where q > 0, r is the time lag and averaging is over the time window T ^> r, usually the 
largest time scale of the system. The function H(q) contains information about averaged 
generalized volatilities at scale r (only q — 1, 2 are used to define the volatility). In particular, 
the Hi exponent indicates persistent (Hi > 1/2) or antipersistent (Hi < 1/2) behavior of 
the trend. For the BRW (brown noise, l// 2 ) one gets H q = 1/2, while for the pink (1//) 
and WN we have H q = 0. For the popular Levy stable and truncated Levy processes with 
parameter a it has been found that H q = q/a for q < a and H q = 1 for q > a |32] . 



The direct computation of the generalized Hurst exponents was performed for return 
times r = 1, 10 and 60 min. as well as for r = 1 day (as for the generalized Renyi exponents 
of the graphs) and for the index q = 0.5, 1.0, 2.0 and 4.0. The value q = cannot be used, 
as q > (see eq. (0)). Computation of H q for the DAX index is displayed in Fig. |9|. The 
linear scaling in the log-log plot is here excellent and the corresponding results are collected 
in Table [TI]. One can clearly see that H q decreases with increasing q, though this effect is 
smaller for longer time lags. Also, for the smaller g-values it decreases with the increasing 
return time r, at least within the range r — 1 4- 30 min., and the process is persistent 
(Hi > 1/2). Very interestingly however, for the larger q- values this tendency gets reversed: 
Hi(r) can be seen to approach the BRW limit from below assuming values significantly 
lower than 1/2 at the high-frequencies. For r ~ 10 min. the exponent Hi approaches the 
BRW value, 1/2. Notice, that at around r c ~ 30 min. the tail index is also stabilized (see 



Sec. II). Similar critical value of r was also reported for S&P500 index in |TTJ and estimated 
t c ~ 20 min. (from (0) it is clear that the structure function S^t) of the index is equivalent 
to the volatility v (r)). Hence, the correlations cease to exist slightly faster for DAX. This 
may indicate that the correlations are weaker in the smaller market. Also, for high frequency 
data the g-dependence is much stronger. When r grows, crossing the critical value r c , H q 
becomes almost independent of q. 

For completeness we have performed the same computation for the DAX returns. The 
quality of the linear scaling in this case is not as excellent as for the index itself but still 
it quite consistently points to H q ~ 0.0 for all the calculated values of q and r, the same 
result as for the WN (in contrast to the index that considerably differs from BRW). This 
difference may be explained by the existence of stronger correlations for the index. 



V. DISCUSSION AND SUMMARY 

The financial processes are governed by a complex dynamics with many degrees of free- 
dom and various additional noise terms due to complicated interactions with the external 
environment. This is a non-stationary evolutionary process with strong correlations, to large 
extent resembling biological evolution, and it cannot be viewed as an equilibrium process 
f33fl . From this perspective it is not very surprising that the formalism that is so far available 
may not yet offer an optimal scheme to consistently comprise all the related effects. The 
present paper constitutes an attempt to provide some further relevant empirical character- 
istics for another important world stock market, the Deutsche Aktienindex, in addition to 
what already is available in the literature. Several such characteristics, not always easy to 
interpret, are identified. 
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In particular, the returns' pdf displays Pareto (fat) tails of the type l/x a . As was 
shown in Sec. II the distribution of the DAX returns has its central part close to the Levy 
distribution with a from 1.2 (minutely returns) up to a ~ 1.7 (daily returns). Within 
the "tail region", i.e. for returns greater than 3 but below 10 4- 20 standard deviations, 
a simple linear regression implies the tail index a ~ 2.4 for minutely returns and up to 
a ~ 3.5 for daily returns. The available data sets seem to be too small to calculate (with 
reasonable certainty level) the tail behavior beyond 10 4 20 standard deviations. We find 
various tail indices for various regions of the distribution, for different return times and for 
different financial indices (see e.g. JTTJ] for comparison). There is no doubt that the power 



like behavior plays important role, but its manifestations are more complex than it was 
originally assumed. It is worth to mention that similar multiscaling of pdf was found in a 
monofractal toy model discussed by Bouchaud et al |24[] . 



From the dynamical point of view more important than pdf are Hurst and Renyi ex- 
ponents of the index and returns' graph. This is because they take into account the time 
dependence of the process. Therefore, the main part of our analysis was devoted to study 
these quantities. For regular (mono-) fractal set the following simple relation between fractal 
exponent of the graph (do) and the Hurst exponent is often satisfied [p], p9| , p0| , p4 | 



do = 1 + E - Hi (8) 

where E is the embedding dimension of the data series (E = 1 for a one dimensional data). In 
particular, this relation is preserved for the WN, for the BRW as well as for the Weierstrass- 



Mandelbrot fractal ]3a] and for any fractional Brownian motion |IJ . This relation is usually 
used when discussing the fractal behavior of financial indices. However, the often cited 
relation (|8|) is not fulfilled in our case. In general, the multifractal formalism for a function, 
instead of the exact spectrum, yields an upper bound of its Holder spectrum [0,0]. Hence, 
in principle it is necessary to calculate the Renyi exponents by direct application of the box 
counting algorithm for the graphs. 

We have used the (2-dimensional) box counting algorithm to calculate the fractal proper- 
ties of the financial graphs. The Renyi exponents with q = 0.5, 1.0, 2.0 and 4.0 for the DAX 
index were found to be around 1.3, slowly decreasing with growing q (about 15% below the 
constant BRW value, 1.5). However, for the DAX returns the corresponding values turn out 
not to be so precisely determined due to a somewhat poorer scaling but still can be localized 
within the range 1.4 4- 1.7, more below the WN value and with stronger g-dependence. The 
decrease of d q with growing q, usually viewed as a sign of multifractality, is bigger here than 
for the index itself, especially for the high-frequency data (1 min. returns). It is worth to no- 
tice that the estimated difference between the fractal exponents of the index and exponents 
of the corresponding differential series (returns) is only about 0.2 4 0.3, while it is exactly 
equal to 0.5 for the BRW and WN. As a potential most visible and significant difference 
between the behavior of the financial index and BRW, this effect demands however further 
more systematic study. Comparing the above results with the calculated Hurst exponents 
one can also see that the relation (H) is violated for the index and seems to be violated even 
more for the returns. In contrast to the Renyi exponents, the generalized Hurst exponents 
for returns are compatible with those of the WN (equal to zero). For the index we obtain 
persistent behavior (Hi = 0.51 4 0.63) which is slightly stronger for short time lags. Here, 
the g-dependence is also stronger than for the d q . 
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One can conclude that the DAX data series, both for the (logarithm of) index and its 
returns, have very complicated self-similar structures that may escape any unique multi- 
fractal description. One likely reason is that the real financial data series are superpositions 
of series with different properties, e.g. for the "regular" and "rush" hours as was suggested 
recently |H-I7[]. To clarify these intriguing effects further empirical and theoretical analysis 
is necessary. 
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TABLES 



TABLE I. Renyi box counting exponents (d q ) for DAX index 



r 


q = 0.5 


q=l 


q = 2 


q = 4 


1 min. 


1.38 


1.36 


1.34 


1.30 


10 min. 


1.36 


1.35 


1.33 


1.29 


60 min. 


1.30 


1.30 


1.29 


1.26 


1 day 


1.33 


1.33 


1.32 


1.29 



TABLE II. Hurst exponents (H q ) for DAX index 



r 


q = 0.5 


q=l 


q = 2 


q = 4 


1 min. 


0.66 


0.63 


0.56 


0.33 


10 min. 


0.52 


0.51 


0.48 


0.36 


60 min. 


0.57 


0.55 


0.51 


0.44 


1 day 


0.52 


0.52 


0.51 


0.46 
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FIGURES 




FIG. 1. Histogram for DAX normalized minutely returns. Normalized Gauss (dotted line) and 
best fit for Levy stable distribution (dashed line, a = 1.16 and 7 = 0.286) are also displayed. 




FIG. 2. Histogram for DAX normalized 10 min. returns. Normalized Gauss (dotted line) and 
best fit for Levy stable distribution (dashed line, a = 1.44 and 7 = 0.327) are also displayed. 
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FIG. 3. Histogram for DAX normalized 1 hour returns. Normalized Gauss (dotted line) and 
best fit for Levy stable distribution (dashed line, a = 1.40 and 7 = 0.354) are also displayed. 
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FIG. 4. Histogram for DAX normalized daily returns. Normalized Gauss (dotted line) and best 
fit for Levy stable distribution (dashed line, a = 1.7 and 7 = 0.385) are also displayed. Lower 
panel: magnification of the tail behavior. 
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og 2 g log 2 g 

FIG. 5. Log-log plot for normalized minutely DAX returns for return times = 1, 10, 60 min. 
and 1 day, respectively for plots (A-D). The triangles represent positive tail while the inverted 
triangles, the negative tail. The dashed line corresponds to the tail index a = 2.4, 2.9, 3.5 and 3.5 
for graphs (A-D). 
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FIG. 6. Recurrence plot for DAX 10 min. returns. Time is on both axis, the embedding 
dimension is set to E = 5 and the normalized neighborhood size is e = 0.2. 
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og 2 N log 2 N 

FIG. 7. Calculating Renyi box counting exponents, d q , for log of DAX. q = 0.5, 1, 2, 4 for 
crosses, triangles, circles, and diamonds, respectively. Plots from A to D are given for return times 
r = 1, 10, 60 min. and 1 day. Solid lines represent least square linear fits. 
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FIG. 8. Calculating Renyi box counting exponents for DAX returns, q = 0.5, 1, 2, 4 for crosses, 
triangles, circles, and diamonds, respectively. Plots from A to D are given for return times r = 1, 
10, 60 min. and 1 day. Solid lines represent least square linear fits. 
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FIG. 9. Calculating generalized Hurst exponents for DAX index, q = 0.5, 1, 2, 4 for crosses, 
triangles, circles, and diamonds, respectively. Plots from A to D are given for return times r = 1, 
10, 60 min. and 1 day. Solid lines represent least square linear fits. 
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